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Preface 

My  main  objective  in  this  study  was  to  develop  a  sub¬ 
sonic  three-degree-of-freedom  store  separation  prediction 
method  which  could  incorporate,  when  available,  experimen¬ 
tally  determined,  free-air,  store  static  stability  data. 

I  am  very  satisfied  with  the  results  I  have  obtained  with 
this  method  and  I  am  hopeful  that  this  longitudinal  analysis 
of  the  separation  problem  will  serve  as  a  precursor  to  a 
six-degree-of -freedom  analysis. 
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assistance  I  received  during  the  course  of  this  study.  I 
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Dr.  Milton  Franke,  for  his  advice  and  guidance.  I  would  also 
like  to  acknowledge  my  other  two  faculty  advisors.  Major  Carl 
Stolberg  and  Captain  James  Karam  f^r  their  constructive 
criticism.  I  owe  a  special  debt  of  gratitude  to  Mr.  Cal  Dyer 
and  Mr.  Jerry  Jenkins,  both  of  AFFDL/FGC,  for  introducing 
me  to  the  store  separation  problem  and  assisting  me  in  the 
computation  and  analysis  of  theoretical  data  from  existing 
methods.  Finally,  I  would  like  to  express  my  deepest  appre¬ 
ciation  to  my  wife,  Ruth,  for  her  patience  and  understanding, 
and  to  my  daughters  Jenny  Lynn  and  Jill  Kathleen  for  the  time 
I  took  from  them  during  these  past  twelve  months. 
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Ab s  tract 

A  subsonic  t hree -d e gre e - of  -  freedom  method  for  predicting 
store  separation  trajectories  is  presented.  The  method  com¬ 
bines  the  calculation  of  interference  loading  on  aircraft 
stores  due  to  F.  Dan  Fernandes  with  loads  due  to  free  air, 
ejector,  and  aerodynamic  damping*  Discussion  of  these  loads 
is  given,  along  with  a  presentation  of  the  three -degree-of - 
freedom  equations  of  motion  and  an  approach  to  their  solution. 
A  computer  program  is  developed  and  used  to  calculate  tra¬ 
jectories  for  the  M-117  bomb  as  carried  on  the  inboard  pylon 
of  the  F-4E.  Comparisons  of  these  trajectory  profiles  with 
wind  tunnel  captive -trajectory  test  results  and  existing 
theory  are  made.  A  user's  guide  and  computer  listing,  ex¬ 
cluding  those  subroutines  due  to  F.  Dan  Fernandes,  are  giver.. 
Recommendations  are  made  for  possible  improvements  to  the 
newly  developed  method. 
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A  SUBSONIC  METHOD  FOR  PREDICTING  THE 
SEPARATION  TRAJECTORIES  OF  AIRCRAFT 
STORES  IN  THE  PITCH  PI.AaE 

I  .  Introduction 

Problem 

Today's  aircraft  are  required  to  carry,  and  subsequently 
to  release,  many  types  of  stores  under  many  different  flight 
conditions.  The  first  0.5  to  1.0  seconds  is  usually  the 
critical  period  in  determining  a  successful  ejection  under 
a  given  flight  condition.  An  analytical  method  of  predicting 
store  separation  trajectories  during  this  critical  time 
period  is  needed  to  serve  as  a  preliminary  design  tool  and 
to  augment  existing  prediction  methods. 

Background 

Store  separation  has  been  a  matter  of  concern  since  the 
first  World  War  I  pilot  attempted  to  throw  a  projectile  from 
his  open -canopied  aircraft  onto  his  enemy  below.  However, 
because  of  the  relatively  slow  speeds  of  early  flight, 
clearance  from  parent  aircraft  posed  no  great  threat  until 
the  advent  of  the  jet.  At  this  time  in  aerial  history, 
pilots  were  required  to  release  externally  carried  stores 
from  aircraft  flying  at  higher  and  higher  speeds  (Ref  12:1,2). 

Prediction  techniques  which  have  been  developed  are 
generally  divided  into  three  categories:  full-scale  flight 
testing,  wind  tunnel,  and  theoretical  methods.  Full-scale 
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flight  testing  has  an  obvious  drawback  in  high  costs,  but, 
more  importantly,  can  be  quite  hazardous  to  the  pilot  and 
his  aircraft.  Wind  tunnel  methods,  including  the  simulated 
trajectory  (Ref  2)  and  grid  techniques  (Ref  1),  have  proved 
to  provide  satisfactory  results.  However,  because  the  only 
Air  Force  tunnel  equipped  for  these  techniques  is  scheduled 
many  months  in  advance  and  because  the  already  expensive 
costs  of  operating  any  wind  tunnel  is  increasing,  these 
methods  cannot  always  be  applied  to  present-day  problems. 

Theoretical  methods  have  only  recently  been  considered 
as  an  alternative  to  the  first  two  approaches.  While  many 
authors  have  developed  analytical  methods  of  calculating 
flow  fields  undep  an  aircraft  (Ref  12:5-19),  Goodwin,  Nielsen 
and  Dillenius  (Ref  7)  and  Goodwin,  Dillenius,  and  Nielsen 
(Ref  8)  have  applied  these  analytical  methods  in  developing 
the  most  comprehensive  technique  for  computing  separation 
trajectories.  Fernandes  (Refs  3  and  4),  has  devised  methods 
for  determining  the  interference  loading  on  an  aircraft 
store  in  the  flow  field  of  a  parent  aircraft.  His  develop¬ 
ment  is  incorporated  in  the  present  work. 

Objective 

The  objective  of  this  study  is  to  develop  a  three- 
degree  -of  -freedom  method  to  predict  the  separation  trajectory 
of  an  aircraft  store  upon  release  from  a  parent  aircraft 
flying  at  subsonic  speeds.  While  previously  developed  tech¬ 
niques  (Refs  7  and  8)  have  been  entirely  analytical  in  nature 
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a  goal  of  this  study  is  to  permit  the  use  of  experimental 
static  stability  data,  when  available,  to  determine  the 
free-air  loads  on  the  store. 

Another  intent  of  this  study  is  to  allow  the  simulation 
and,  therefore,  the  integrated  effects,  of  a  two- point 
ejector  system  in  calculating  t h 2  trajectory  of  an  aircraft 
store.  This  again  contrasts  previously  developed  techniques 
which  do  not  simulate  the  ejector,  but  account  for  it  in 
initial  conditions  which  must  be  calculated  before  trajec¬ 
tory  calculations  can  begin. 

A  shortcoming  of  previously  developed  methods  is  the 
large  computer  time  necessary  to  simulate  the  separation 
trajectories  of  pircraft  stores.  A  third  purpose  for  under¬ 
taking  this  study  is  to  write  a  FORTRAN  computer  program 
which  will  require  reasonable  computer  execution  time  to 
simulate  the  trajectories  of  ejected  stores. 

Finally,  it  is  expected  that  the  t hr ee -d egree -o f -freedom 
store  separation  prediction  method  will  serve  as  a  precursor 
to  a  si x -degree -of - fr eed om  method. 

Approach 

The  method  used  to  calculate  separation  trajectories 
is  straightforward.  Interference  normal  forces  and  pitching 
moment  coefficients  on  the  ejected  store  are  determined 
using  portions  of  Fernandes*  subsonic  interference  loading 
program  (Ref  3).  Predetermined  free-air  normal  force  and 
pitching  moment  coefficients,  obtained  either  by  experiment 
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or  any  suitable  theory,  are  then  added  to  the  interference 
loads.  This  sum,  along  with  a  predetermined  axial  force 
coefficient  which  remains  constant  throughout  the  trajectory, 
constitutes  the  total  aerodynamic  loading  on  the  store. 

Next,  normal  force  and  pitching  moment  coefficients  due  to 
the  store  ejector  system,  when  operating,  are  determined. 
Finally  the  pitching  moment  coefficient  due  to  damping  is 
calculated  using  a  predetermined,  from  experiment  or  theory, 
pitch  damping  coefficient.  All  force  and  moment  coefficients 
are  then  summed  and  a  new  position  is  determined  by  inte¬ 
grating  the  t hree -de gre e -o f - freed om  equations  of  motion. 

The  process  is  repeated  until  clearance  from  the  parent  air¬ 
craft  is  assured. 
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II.  Discussion  of  Hoads  o n  t he  Store 

Interference  Loads 

Interference  pitching  moment  and  normal  force  coeffi¬ 
cients  are  calculated  using  portions  of  the  F.  Dan  Fernandes 
subsonic  interference  loading  program  (Ref  3).  The  calcula¬ 
tion  of  the  interference  flow  field  by  this  method  entails 
a  linear  theory  of  source,  vortex,  and  doublet  distributions. 
The  aircraft  wing,  pylon,  inlet  and  fuselage  nose  are  modeled 
using  various  combinations  of  these  distributions,  the 
strengths  of  which  are  calculated  by  satisfying  boundary 
conditions  (no  flow  normal  to  the  body  surface)  at  various 
control  points  on  the  body.  Disturbance  velocities  are  then 
calculated  over  the  length  of  the  store  as  functions  of 
these  strengths  and  distances  to  the  field  point  under  con¬ 
sideration.  From  these  disturbance  velocities  an  interfer¬ 
ence  angle-of-attack  field  can  be  determined.  Interference 
forces  and  moments  are  then  calculated  by  integrating  the 
effects  of  this  variable  interference  angle-of-attack  field 
over  the  length  of  the  store  using  predetermined  (by  experi¬ 
ment  or  theory  provided  by  Fernandes)  free-air  body  loading 
coefficients  per  unit  angle-of-attack  per  unit  length. 
Concurrently,  the  interference  static  pressure  field  over 
the  length  of  the  store  is  calculated  and  then  integrated 
to  yield  the  loading  due  to  buoyancy  on  the  store  (Ref  3:5, 
E-l).  Compressibility  effects  are  included  by  applying  the 
Pr and t 1 -G 1 auert  rule  (Ref  3:A-1).  Details  on  the  interference 
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loaJing  method  are  given  by  Fernandes  (Ref  3). 

Although  Fernandes'  program  calculates  five  interfer¬ 
ence  force  and  moment  coefficients  (normal  force,  sid' 
force,  pitching  moment,  rolling  moment,  and  yawing  moment 
coefficients),  the  present  method  uses  only  two  (pitching 
moment  and  normal  force  coefficients)  because  it  is  a 
three -degree -of- freedom  analysis  only.  Interference  axial 
force  coefficient  is  not  calculated  by  the  Fernandes  method 
and  so  this  interference  force  is  not  accounted  for  in  the 
present  method  . 

Free -Air  Load s 

The  calculation  of  total  static  aerodynamic  loads  on 
the  store  is  completed  with  the  addition  of  free-air  loads 
to  those  interference  loads  determined  by  the  Fernandes 
method.  Free-air  normal  force  and  pitching  moment  coeffi¬ 
cients  are  not  calculated  by  the  present  method  and,  there¬ 
fore,  must  be  obtained  before  using  it,  usually  as  functions 
of  store  angle  of  attack  and  store  Much  number.  The  free- 
air  axial  force  coefficient  must  also  be  obtained  before 
using  the  present  method  and  is  r  turned  constant  throughout 
the  trajectory  simulation.  Data  may  be  obtained  from  experi¬ 
mental  sources  (as  was  done  for  the  sample  case  considered 
in  this  report--see  Chapter  V)  or  from  any  suitable  analyti¬ 
cal  method,  such  as  USAF  DATCOM  (Ref  10,  Sec  4). 
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I:  j  c c t  o r  Load s 

Calculation  of  ejector  normal  forces  and  pitching 
moments  is  quite  similar  to  the  method  used  by  Christopher 
and  Carlcton  (Ref  2:32-33).  Modifications  were  necessary 
to  allow  ejector  force  curves  either  as  a  function  of  ejector 
foot  distance  or  as  a  function  of  time.  In  addition  calcula¬ 
tion  of  ejector  normal  force  and  pitching  moment  coefficients 
are  made  to  permit  comparisons  of  those  ejector  coefficients 
with  corresponding  interference  and  free-air  coefficients. 

The  present  method  requires  the  ejector  force-distance 
or  force-time  curves  to  be  in  the  form  of  a  fifth  degree 
polynomial,  the  same  form  used  by  Christopher  and  Carleton 
(Ref  2:32).  Two  such  curves  are  allowed  to  accommodate 
a  two-point  ejector  system.  In  addition,  the  distance  away 
from  the  center -of -grav it y  of  the  store  at  which  each  ejector 
is  assumed  to  act  must  be  known  to  calculate  ejector  moments 
on  the  store. 

Pitching  Moment  Due  to  Damping 

To  account  for  the  change  in  pitching  moment  due  to 
aerodynamic  damping,  a  pitch  damping  coefficient,  required 
by  the  present  method,  must  be  known.  According  to 
Christopher  and  Carleton,  the  consideration  of  pitching 
moment  due  to  aerodynamic  damping  is  usually  a  second  order 
effect,  so  an/  reasonable  estimate  for  pitch  damping  coeffi¬ 
cient,  C  ,  should  satisfy  the  requirements  of  the  present 
method.  Should  the  store  exhibit  large  oscillatory  motions. 
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however,  this  effect  could  become  significant  in  which  case 
consideration  should  be  given  to  obtaining  experimental 
values  for  C^,  or  running  trajectory  sensitivity  studies 
varying  possible  values  for  Cm^  in  order  to  determine  its 
significance  on  the  trajectory  of  the  store  under  study 
(Ref  2:9) . 

The  actual  calculation  of  pitching  moment  coefficient 
due  to  aerodynamic  damping  is  a  simple  one  once 


C  .  dC/d(qb/2  V, 


J  ,  is  known  : 


=  C 


DAMP 


mq 


qb/2  V 


store 


where 


q  =  store  pitch  rate 
b  =  store  reference  length 
'store  =  abs°lute  store  total  velocity, 


Total  Load  s  on  the  Store 

The  total  normal  force  load  on  the  store  is  given  by 
the  sum  of  the  loads  due  to  interference,  free-air,  and 
ejector  systems.  The  total  axial  force  load  on  the  store 
is  computed  using  the  predetermined  axial  force  coefficient 
The  total  pitching  moment  on  the  store  is  simply  the  sum  of 
the  moments  due  to  interference,  free-air,  ejector  system, 
and  aerodynamic  damping.  These  loads  are  illustrated  in 
Fig.  1,  in  the  coordinate  system  fixed  in  the  ejected  store 


(1) 
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V/HERE 

fAERO,N  =  sun  of  interference  and  free-air  aerodynamic 
forces  in  the  normal  direction 

^AEROjA  =  free-air  aerodynamic  force  in  the  axial  direction 

Fzl  =  force  due  to  ejector  one 

pz2  =  force  due  to  ejector  twt 

mAERO  =  pitching  moment  due  to  interference  and  free- 
air  aerodynamic  loads 

NDAMP  =  pitching  moment  due  to  damping 
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III.  Integration  o  f  Forces  nml  Moments 

i 

liquations  o  f  Motion 

The  form  of  the  thrce-degrce-of-freedoin  equations  of 
notion  presented  in  this  chapter  was  taken  from  Goodwin, 

Nielsen,  and  Dillenius  (Ref  7:78-84).  In  place  of  the  Goodwin 
method  for  calculating  forces  and  moments  on  the  store,  how¬ 
ever,  the  calculation  of  forces  and  moments  by  the  present 
method,  including  free-air  loading,  interference  loading, 
aerodynamic  damping,  and  ejector  loading,  is  used. 

The  inertial  f  t  ame  (£  ,  p )  for  this  treatnent  is  fixed 
with  the  aircraft,  hence  the  flight  conditions  must  be  non¬ 
accelerating.  The  aircraft,  therefore,  is  assumed  to  fly  at 
a  constant  flight  path  angle,  y_,  free  stream  velocity,  V  , 
and  angle  of  attack,  a„ .  Figure  2  illustrates  this  inertial 

D 

coordinate  system  and  the  initial  orientatirn  of  the  store 
in  it . 

The  longitudinal  equations  of  motion  in  this  coordinate 
system  are: 


mj  =  F- 

(2) 

s  <; 

m  n  =  F 

(3) 

s  n 

m  k20  =  Mn 
s  0 

(4) 

Now  the  forces  on  the  store  in  the  ^-direction  and 
q - di rect i on ,  and  F^,  can  be  found  by  properly  resolving 
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forces  along  the  store  axis  system  into  the  inertial  frame. 
The  force  normal  to  the  store  is  simply 


F,,  =  F .  .  +  F  r  +  F  .  +  F  0 
h  int  fa  zl  z2 


(5) 


The  axial  force  on  the  store  is  given  by 


FA  “  2  ^“^store^ref 

The  resolved  forces  are,  then 

F^  =  Fn  sin  0  +  FA  cos  (CVYB-YS)  +  msg  sin  (ag+Yg) 
Fn  =  FN  C0S  6  +  FA  Sin  (aB  +  YB'ys:i  -  mS2  cos  (aB  +  YB) 


(6) 

(7) 

(8) 


where  Yg  is  The  flight  path  angle  of  the  store,  which  is  dif¬ 
ferent  from  yb  if  wind  tunnel  captive  trajectory  tests  are  to 
be  simulated  (Ref  7:84-85). 

Noting  that  CN  =  FN/(~  PcoVstoreSref 3  and  substitutin8 
Eqs  (7)  and  (8)  into  Eqs  (2)  and  (3),  the  first  two  equations 
of  motion  become 


£ 


T  Poo^store  ^re f'  "‘s  '  lv>N 


:/»0  lCM  sin 


0  +  C.  cos 
A 


^(2B  +  YB"YS^  ^ 


+  g  sin (ag+Yg) 


(9) 


=  I  P»Vstora(Sref/ms)[CN  005  9  *  CA  sin  'VW1 


-  g  cos(aB+YB) 


(10) 


Substitution  of  C  =  Mn/(4  p  S  rb)  into  Eq  (9) 

m  0  2  °°  store  ref  n 

yields  the  third  longitudinal  equation  of  motion 
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^v;tore(Srefb^sk~)C, 


(11) 


where  C  is  the  sun  of  the  interference,  free- air,  eiector 
rn 

one,  eiector  two,  and  pitch  d nr. ring  r.orcnt  coefficients. 
Equations  (9),  (10),  and  (11)  are,  except  for  slight  changes 

in  notation,  exactly  those  equations  derived  by  Goodwin, 
Nielsen,  and  Dillenius  (Ref  7:79-80). 


Integra tion  Techniques 

Equations  (9),  (10),  and  (11)  nust  be  integrated  in 
order  to  calculate  new  store  positions  as  well  as  angular 
and  linear  velocities  and  accelerations  of  the  store  at  each 
tine  of  interest.  The  nethod  used  for  integration  is  a 
standard  fourth-order  predictor-corrector  technique,  utilizing 
a  Runge-Kutta  schene  to  calculate  intermediate  steps,  and 
was  taken  from  Goodwin,  et  3l.  The  method  is  introduced  in 
their  t’nre  e  -  de  g  re  e  -  o  f  -  f  reedon  report  (Ref  6:26)  and  explained 
in  detail  in  their  six-degree-of-freedon  report  (Ref  9:139-142, 
219-221)  . 

Actual  integration  of  the  equations  of  notion  nay  be 
started  at  any  desired  time,  provided  the  correct  initial 
conditions,  store  position  and  linear  and  angular  velocities, 
are  known.  The  forces  and  moments  are  calculated  using 
methods  outlined  in  Chapter  II,  then  a  new  position  is  found 
by  integrating  Eqs  (9),  (10),  and  (11)  into  which  have  been 
substituted  the  calculated  values  for  C  and  C„.  The  process 
is  then  repeated  until  clearance  is  assured. 
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I V .  Computer  Program 

Computer  Memory  and  execution  T i me  Requi renents 

A  FORTRAN  computer  program  was  developed  for  the  present 
method.  It  is  operational  on  the  CDC  6600  computer  with  the 
SCOPE  3.4  operating  system  and  library  tape  and  requires  about 
60000  octal  storage  registers  for  loading.  Execution  time 
varies  greatly,  depending  on  how  much  detail  is  taken  in 
describing  the  geometry  of  the  parent  aircraft,  how  many  store 
sections  are  considered  for  computing  interference  loading, 
the  time  step  size  chosen  for  integration  of  the  equations 
of  motion,  and  the  computational  mode  desired  for  this 
integration.  Execution  time  for  the  trajectory  simulations 
of  the  M-117  bomb  from  an  F-4E  are  listed  in  Chapter  V  and 
discussion  of  computational  modes  in  Appendix  A. 

Overview  of  Computer  Program 

The  bulk  of  the  computer  program  consists  of  twenty- six 
subroutines  taken  from  Fernandes  (Ref  5),  all  necessary  in 
computing  the  interference  loading  on  the  aircraft  stores. 

In  addition,  the  predictor-corrector  subroutine  necessary 
for  integration  of  the  equations  of  motion,  subroutine  ADAMS, 
was  taken  from  Goodwin  (Ref  6). 

The  user  of  the  computer  program  must  supply  his  own 
subroutine,  FREAIR,  whose  purpose  is  to  calculate  the  free- 
air  pitching  moment  and  normal  force  coefficients  on  the 
store,  given  the  store  angle  of  attack  and  Mach  number. 
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The  main  program  whose  primary  purpose  is  to  organize 
calculations  by  calling  necessary  subroutines,  calls  two 
other  subprograms,  TRREAD  and  TRAJEC.  The  purpose  of  sub¬ 
routine  TRREAD  is  to  read  in  values  of  store  geometry,  store 
mass  characteristics,  ejector  characteristics,  pitch  damping 
coefficient,  axial  force  coefficient,  and  aircraft  flight 
conditions  necessary  for  separation  calculations.  This  sub¬ 
routine  also  initiates  values  of  the  dependent  variables, 
store  position  and  pitch  angle,  and  linear  and  angular  veloc¬ 
ities  required  for  integration  of  the  equations  of  motion. 

The  purpose  of  subroutine  TRAJEC  is  to  accept  store 
position  and  interference  loads,  calculate  free-air  loads, 
ejector  loads,  and  pitch  damping,  sum  all  loads  and  output 
them.  It  then  integrates  the  equations  of  motion  and  returns 
a  new  store  position  and  orientation  to  the  main  program. 

A  user's  guide  for  the  computer  program  is  presented  in 
Appendix  B  and  a  listing  of  the  program,  with  the  exception 
of  those  26  subroutines  taken  from  Fernandes  (Ref  5),  is  pro¬ 
vided  in  Appendix  C.  A  few  minor  changes  to  those  Fernandes 
subroutines,  necessary  to  pass  interference  coefficients  to 
the  main  program  and  to  avoid  extraneous  output,  are  also 
listed  in  Appendix  C.  The  computer  codes  for  the  Fernandes 
subroutines  themselves  are  available  through  COSMIC  (Computer 
Software  Management  and  Information  Center).  Requests  should 
be  directed  to:  COSMIC,  University  of  Georgia  30601. 
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Parameter  Name 

Parameter  Value 

ms  >  slugs 

25.45 

XCG,  feet 

2.333 

X^,  feet 

-0.3117 

Z  ,  feet 

Lmax 

0.275 

Sref’  Sq  ft 

1  .  395 

b,  ft 

1 .333 

n 

> 

0.1 

IyY»  slugs-sq  ft 

5  0.0 

C  ,  per  rad 
mq  r 

-70.0 

n  for  a  sample  calculation  be- 
ee-air  stability  characteristics 
d,  under  certain  flight 
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Sketch  of  the  M-117  All-Purpose  Bomb 
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conditions  and  ejector  loads,  the  M-117  exhibited  small 
lateral  (yaw  angle  and  y-direction)  excursions. 

Experimental  trajectories  were  taken  from  an  A II DC  wind 
tunnel  report  whose  author  investigated  the  effects  of  wing 
leading-edge  slats  on  the  separation  characteristics  of  vari¬ 
ous  stores  as  carried  on  the  F-4E’  (Ref  11).  That  same  wind 
tunnel  investigation  modeled  separation  trajectories  from  the 
F-4E  aircraft  without  leading  edge  modifications  and  it  is 
those  baseline  trajectories  to  which  the  th'oretical  results 
are  compared. 

The  ejector  force-distance  curve  used  in  the  present 
method  was  taken  from  the  same  AF.DC  wind  tunnel  report  and 
is  shown  in  Fig.  4. 

Comparison  of  Store  Trajectories 

Figures  5,  6,  7  and  8  exhibit  X  and  Z  coordinates,  rela¬ 
tive  to  carriage  position,  of  the  center  of  gravity  of  the 
store,  as  well  as  the  absolute  pitch  angle  of  the  store  as 
seen  in  the  fuselage  wind-axis  system.  Results  from  experi¬ 
ment  and  both  the  present  method  as  well  as  the  method  due 
to  Dillenius,  et  al.  (Ref  8)  are  shown. 

As  illustrated  in  Fig.  5  both  methods  predicted  the  ex¬ 
perimental  results  remarkably  well  for  the  low  speed,  moderate 
angle-of- attack  esse  (Mach  =  0.332,  a  =  7.40),  with  the 
present  method  doing  slightly  better  in  position  prediction 
and  the  previously  developed  method  predicting  pitch  angle 
excursions  more  accurately.  In  Fig.  6,  however,  while  linear 
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Fig.  8.  Simulated  Separation  Characteristics  of 
the  M-117  Bomb. 
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distance  vs.  time  curves  are  a  5;  a  in  estimated  quite  well,  the 
pitch  angle  calculation  is  not  as  close  -to  experiment  as  in 
the  previous  case.  It  is  difficult  to  conjecture  which 
method  predicts  pitch  angle  better,  since  both  indicate  the 
damped  sinusoidal  pitch  angle  profile  common  to  stable 
stores,  but  the  trajectory  due  to'  the  present  method  has  a 
shorter  period  of  pitch  oscillation  than  does  the  trajectory 
due  to  the  Goodwin  method  (Ref  8). 

For  the  high  speed  case  (Mach  =  0.829),  F’gs.  7  and  8 
show  interesting  results.  While  the  present  method  again 
predicts  linear  travel  quite  well  and  pitching  angle  reason¬ 
ably  well,  the  method  obtained  from  Ref  8  shows  a  definite 
deviation  from  experiment  slightly  after  ejector  time  cut-off 
(about  0.08  seconds).  A  Mach  number  equal  to  0.829  results, 
more  than  likely,  in  a  supercritical  flow-field  situation, 
thus  violating  a  basic  assumption  of  the  method  (Ref  8),  that 
conditions  modeled  are  to  be  at  subcritical  speeds  only, 
thereby  rendering  the  results  calculated  at  the  high  speeds 
invalid. 

The  application  of  the  present  method  to  relatively  high 
Mach  numbers  (0.829)  might  also  be  construed  as  a  violation 
of  the  subcritical  speed  restrretion  placed  on  the  Fernandes 
method.  While  technically  tnis  is  true,  there  is  no  such 
restriction  on  the  other  sources  of  store  loading,  so  the 
present  method  can  be  pushed  past  the  normal  critical  speed 
(about  Mach  =  0.7  to  0.8)  with  relative  safety  provided  the 
free-air  loading  is  known  with  some  confidence. 
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Reasons  for  the  excellent  agreement  of  the  present 
method  with  experiment  are  many,  one  of  which  is  certainly 
the  acceptable  results  the  Fernandes  method  for  calculating 
interference  forces  and  moments  seems  to  produce  (Ref  3:13-21). 
Good  results  for  total  aerodynamic  loading  is  assured  by  the 
fact  that  the  experimental  free-air  data  was  available  to  add 
to  the  analytically  predicted  interference  loads. 

Store  loading  due  to  ejector  forces  and  due  to  aero¬ 
dynamic  damping  are  modeled  in  the  present  work  exactly  as 
they  are  in  the  AH DC  tunnel  tests  and  as  a  result  should 
introduce  no  error  into  calculations  of  store  trajectories. 

In  the  case  of  the  ejector  forces  at  low  speeds,  this  efvect 
can  be  quite  significant  because,  when  operating,  the  ejector 
forces  tend  to  overwhelm  the  aerodynamic  forces.  Figure  9 
illustrates  this  point  for  the  low  speed  (Mach  =  0.332)  case. 
At  higher  aircraft  speeds,  however,  the  aerodynamic  forces 
are  greater  and  do  not,  in  general,  dominate  the  ejector 
forces,  as  indicated  in  Fig.  10  for  the  Mach  =  0.829  case. 

A  final  reason  for  satisfactory  results  is  the  fact  that 
the  M-117  is  a  very  dense  store.  Dense  stores  have  higher 
inert i al  -  to - aerodynami c -  1 o ad s  ratios  than  do  lighter  weight 
stores  of  the  same  general  shape  and  so  are  not  "blown 
around"  as  easily  as  their  lighter  weight  counterparts. 

Comparison  o f  Computer  Execution  Time 

The  execution  times  required  by  the  computer  program 
utilizing  the  present  method  can  be  directly  compared  to  that 
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Fig.  9,  Comparison  of  Loads  on  the  M-117  Bomb  for 
the  Low  Speed  Case. 
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of  the  conputer  program  written  by  Goodwin  and  Dillenius 
(Ref  9).  Two  different  execution  times  are  listed  in  Table  II 
for  the  present  method,  illustrating  the  two  computational 
modes  available  in  its  corresponding  computer  program.  These 

Table  II 


Comparison  o_f  CPC  660  0  Computer  P x e c u t i o n  Time 


Computer  Execution  Time  (secs) 

Present 

Method 

Ref  8 
Method 

Case 

Mach 

Alpha 

KSTABLE=0 

KSTAB LE  = 1 

1 

0.332 

7.4 

98 

30 

205 

2 

0.540 

11.0 

98 

67 

204 

3 

0.829 

0.4 

142 

46 

•  209 

4 

0.829 

3.6 

142 

43 

210 

computational  modes,  KSTABLE=0  and  KSTABLE=1,  yield  essentially 
the  same  store  trajectories  and  are  explained  in  detail  in 
Appendix  A.  Examination  of  Table  II  reveals  that  the  present 
method  yields  some  dramatic  savings  in  computer  execution 
time,  while  losing  no  noticeable  accuracy  in  trajectory  simu¬ 
lation  at  low  and  moderate  speeds,  and  actually  extending  the 
speed  regimes  for  reasonably  accurate  trajectory  simulation 
to  higher  subsonic  speeds  than  those  allowed  by  the  Goodwin 
method  (Ref  8) . 

Both  computer  programs  are  operational  on  the  CDC  6600 
computer  with  the  SCOPE  3.4  operating  system  and  library 
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tape.  The  program  due  to  Goodwin  (Ref  9),  which  was  modified 
to  restrict  movement  to  the  pitch  plane  only,  requires  about 
114000  octal  storage  registers  for  loading.  The  program 
written  using  the  present  method  requires  about  60000  octal 
store  registers  for  loading. 
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VI.  Conclusions  and  Recommendations 


Conclusions 

Results  from  the  present  prediction  method  compared 
well  with  experimental  results  for  store  trajectories  which 
exhibit  minimal  lateral  characteristics.  Although  previously 
developed  trajectory  prediction  techniques  (Refs  7  and  8) 
and  the  method  used  in  the  present  work  to  calculate  inter¬ 
ference  loads  are  limited  to  the  subcritical  speed  regime, 
the  present  method  may  be  used  with  some  confidence  past  the 
subcritical  cut-off  because  experimental  data  used  to  calcu¬ 
late  free-air  loads  on  the  store  is  usually  available  for 
supercritical  speeds.  In  addition,  the  method  used  to  calcu¬ 
late  the  loads  due  to  the  ejector  system  is  not  limited  to 
any  aircraft  speed  regime,  so  the  method  used  to  calculate 
ejector  loads  may  also  be  used  supercri t i ca 1 ly . 

A  substantial  improvement  in  computer  execution  time 
was  realized  with  the  present  method  over  the  method  due  to 
Goodwin,  Dillenius,  and  Nielsen  (Ref  8).  This  savings  was 
even  greater  for  the  KSTABL0=1  mode,  discussed  in  Appendix  A, 
although  all  store  trajectory  simulations  may  not  be  suitable 
to  allow  the  use  of  that  time-saving  feature. 

Regardless  of  how  well  the  present  method  predicts  the 
t hree -de gree - o f - freedom  trajectory  of  an  aircraft  store,  its 
engineering  practicality  is  limited  because  it  is  a  three, 
not  a  six,  degre e- o f - freedom  analysis.  Should  the  store 
trajectory  under  investigation  not  be  suited  to  a  pitch-plane 
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analysis,  or  should  experimental  store  static  stability  data 
not  be  available,  the  six-degree-of-freedom  analysis  due  to 
Goodwin,  Dillenius,  and  Nielsen  (Ref  8)  should  be  used. 

Recommendations 

As  alluded  to  above,  the  first  obvious  extension  of  the 
present  method  would  be  one  to  accommodate  the  other  three 
degrees  of  freedom  (side  force,  yawing  moment,  and  rolling 
moment).  While  this  sounds  straightforward,  consideration 
of  all  six  degrees  of  freedom  compounds  problems  immensely 
and  should  not  be  viewed  lightly.  Goodwin,  Dillenius,  and 
Nielsen  (Ref  8)  have  completed  this  task  using  their  theory 
for  calculating  store  loading,  so  much  of  their  work  could 
again  be  utilized.  However,  as  yet,  the  Fernandes  method 
does  not  allow  consideration  of  a  yawed  store  so  consider¬ 
able  modification  to  his  method  must  be  completed. 

In  addition,  improvement  to  the  Fernandes  method  itself 
could  be  attempted.  Of  particular  interest  would  be  to  allow 
some  consideration  of  wing-body  interference,  addition  of  an 
arbitrarily  shaped  representation  of  the  fuselage,  and  the 
addition  of  other  bodies  to  account  for  store-to-store  inter¬ 
ference.  A  method  similar  to  Goodwin,  et  al .  (Ref  7:20-22) 
would  be  the  simplest  to  employ  in  order  to  get  a  first  order 
effect  for  wing-fuselage  and  store-to-store  interference. 

In  that  method,  an  induced  camber  on  the  wing  is  calculated 
to  cancel  the  effect  of  the  fuselage  changing  the  boundary 
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conditions  on  the  wing  when  the  two  independent  solutions, 
wing  and  hotly,  arc  superimposed. 

Thirdly,  the  rack  could  he  modeled,  possibly  as  a  dis¬ 
tribution  of  axi symmetric  sources  (Ref  8:20)  or,  as  suggested 
by  Goodwin,  et  al.  (Ref  8:20),  a  combination  distribution  of 
axisymmetric  sources  to  account  for  rack  thickness  and  a 
small  system  of  vortices  to  account  for  the  short  wing-like 
stubs  which  normally  protrude  from  most  bomb  racks. 

Another  area  which  typically  gets  little  consideration 
is  the  elastic  effect  on  the  rack  and,  consequently,  the 
net  force  felt  by  the  store  duo  to  the  extremely  high  ejector 
forces  common  to  most  store-ejector  systems. 

Finally,  consideration  should  be  made  to  the  development 
of  a  supersonic  three  or  six-degree -of- freedom  trajectory 
program.  Fernandes  (Ref  4)  has  written  a  companion  program 
to  his  subsonic  interference  loading  program  which  will 
compute  interference  loads  on  a  store  in  the  flow  field  of 
an  aircraft  flying  at  supersonic  speeds.  Ilis  supersonic 
method  could  be  incorporated  into  a  new  technique  in  a  manner 
quite  similar  to  that  developed  here. 
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Appendi x  A 

C o m p » t  n  t  L o n  a  1  Mode s  A v a  i  1  a h  I  c  i n  the 
Conn uter  P  r op. ram 

A  standard  f our t h - o rde r  p r c d i c t or  -  cor rec t o r  technique 
wac  used  to  integrate  the  equations  of  motion.  At  the  onset, 
and  at  each  subsequent  change  in  step  size  of  the  integration 
routine,  Runge-Kutta  calculations  are  made  to  determine 
intermediate  values  of  the  dependent  variables.  Large  cojd- 
puter  execution  times  were  introduced  because  at  each  major 
and  each  intermediate  time  step  all  forces  and  moments  must 
be  calculated  according  to  the  methods  outlined  in  Chapter 
II.  Calculation  of  interference  forces  and  moments,  as  might 
be  expected,  required  the  most  amount  of  time,  on  the  order 
of  two  seconds  per  store  position  for  a  typical  f i ghter- bomber . 
To  alleviate  this  problem  a  parameter,  KSTABLE,  may  be  input 
with  the  value  of  1.  This  will  suppress  the  calculation  of 
interference  loadings  at  each  intermediate  time  step,  thereby 
reducing  computational  time  significantly.  The  resultant 
savings  in  computational  time  for  the  trajectory  simulation 
examined  in  this  study  is  listed  in  Table  II  of  Chapter  V. 

Some  care  should  be  exercised  in  choosing  which  computa¬ 
tional  mode,  the  KSTABLE=0  mode  or  the  KSTABLE=1  mode,  is 
used.  In  the  example  of  the  M-117  bomb  trajectory  simula¬ 
tion,  interference  loads  did  not  change  much  for  intermediate 
time  steps.  This  was  probably  due  to  the  high  weight  and 
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good  stability  of  the  bomb  causing  small  pitch  oscillations 
at  these  intermediate  time  steps.  This  may  not  be  the  case, 
however,  for  unstable  or  lighter  weight  stores,  so  trajectory 
sensitivity  studies  should  be  run,  using  both  computational 
modes,  at  extreme  Mach  numbers  and  angles  of  attack.  A 
comparison  of  trajectories  calculated  from  each  mode  should 
enable  the  determination  of  the  suitability  of  the  time-saving 
KST AB Lb  =  1  mode  for  the  ai rcraft/store  under  consideration. 
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Appendix  B 

User's  Guide  to  the  Conputcr  Program 

Most  of  the  inputs  to  the  conputcr  program  describe 
the  aircraft  geometry,  store  loading,  coefficients,  and 
parameters  to  determine  the  number  of  flow  singularities 
used  to  represent  the  interference  flow  field  as  seen  by 
the  store.  These  items  are  input  exactly  as  described  by 
Fernandes  (Ref  3:26-47)  and  will  not  be  repeated  here.  The 
only  change  from  the  original  Fernandes  inputs  is  his  control 
parameter  for  run  stacking,  IGO.  This  final  Fernandes  input 
is  eliminated  altogether  for  the  present  computer  p'ogram. 

The  remaining  inputs  required  by  the  present  computer 
program  will  now  be  presented. 

CARD  1  (after  Fernandes  inputs)  consists  of  four 
control  parameters  input  in  format  415. 

KSTABLE  is  the  computational  mode  parameter  discussed 
in  Appendix  A.  KSTABLE=0  causes  calculation 
of  interference  forces  at  all  time  steps; 
KSTABLE=1  causes  calculation  of  interference 
forces  at  major  time  steps  only. 

NEJECT  is  the  ejector  curve  control  parameter. 

NJECT=0  implies  the  ejector  curve  to  be  input 
is  For ce - di s t ance  ;  NEJECT=1  implies  Force¬ 


time. 


CiAI.  /  Mi  /  7  4  S  -  4 


r 

NGAM  is  the  index  controlling  flight  path  angle. 

i\’GAM  =  0  causes  the  trajectory  simulation  to  be 
free-air;  N G A M  = 1  causes  a  wind  tunnel  c a p t i v e - 
store  simulation. 

N'CNMAX  is  the  index  indicating  the  maximum  number 

of  major  time  steps  allowed  during  integration. 

CARD  2  contains  foui  store  parameters  input  in 
4F10.5  format. 

SMASS  is  store  mass.  (slugs) 

RGYRAY  is  store  radius  of  gyration  about  y-axis. 

( feet) 

XCG  is  the  store  center  of  gravity  measured  behind 

store  nose.  (feet) 

SRMAX  is  store  maximum  radius.  (feet) 

CARD  3  contains  these  four  parameters  input  in 
4E12.4  format. 

VI NF  is  aircraft  flight  velocity.  (feet  per  second) 

GAMF  is  aircraft  flight  path  angle.  (degrees) 

RHO  is  air  density  at  simulated  altitude.  (slugs 

per  cubic  foot) 

G  is  acceleration  due  to  gravity.  (feet  per 

second  per  second) 

CARD  4  requires  these  four  inputs  in  format  4E12.4. 

CA  is  the  store  axial  force  coefficient,  assumed 

constant  throughout  the  simulation. 


A 
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is  the  store  initial  translational  motion. 


Direction  is  normal  to  store  longitudinal  axis, 
(feet  per  second) 

VAR(6)  is  the  store  initial  pitching  velocity  about 
y-axis.  Positive  is  nose  up.  (radians  per 
second). 

Note  that  vZHRO  and  VAR(6)  should  be  input  as  zero  unless 
ejector  simulation  is  not  desired. 

CARD  5  contains  a  one - d i mens i onal  array  of  order 
six  whose  input  is  6E12.4. 

C(l),  C(2),  ...  C(6)  are  the  coefficients,  low-to-high 

order,  of  the  fifth  order  polynomial  curve 
fit  of  the  ejector  Force-time  or  Force-distance 
curves  of  ejector  one. 

CARD  6  also  contains  a  one-dimensional  array  of 
order  six  whose  input  is  6 F.  12.4. 

D  ( 1 )  ,  D(2),  ...  D ( 6 )  are  the  coefficients,  low-to-high 
order,  of  the  fifth  order  polynomial  curve 
fit  of  the  ejector  Force-time  or  Force-distance 
curves  of  ejector  two. 

CARD  7  consists  of  six  parameters  input  in  format 

6F10.5. 

CMQ  is  the  pitch  damping  coefficient.  (per  radian) 

XL1  is  the  ejector  one  piston  location  relative 

to  the  store  center  of  gravity,  positive 
forward  of  store  center  of  gravity.  (feet) 
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XL2  is  the  ejector  two  piston  location  relative 

to  the  store  center  of  gravity,  positive 
forward  of  store  center  of  gravity.  (feet) 

EJE.ND1  is  the  ejector  one  cut-off  argument.  (feet 
or  seconds,  depending  on  NEJECT) 
nJEN'D2  is  the  ejector  two  cut-off  argument.  (feet 
or  seconds,  depending  on  NEJECT) 

EJANGL  is  the  angle  from  the  vertical  in  the  aircraft 
y-z  plane,  at  which  each  ejector  is  assumed 
to  act.  (degrees) 

Note  that  in  this  pitch-plane  analysis  only  that  component 
of  the  ejector  force  curve  in  the  aircraft  z-dircction  is 
considered . 

CARD  8  contains  four  time  parameters  input  in 
format  4E12.4. 

DTI  ME  is  the  initial  integration  interval.  (seconds) 

TIMEI  is  the  initial  time  of  trajectory.  (seconds) 

TIMEF  is  the  final  time  of  trajectory.  (seconds) 

DTIME2  is  the  integration  interval  after  both  ejectors 

have  stopped.  (seconds) 

If  no  ejector  is  present  input  DTIME2  equal  to  DTIME. 

If  at  least  one  ejector  is  present,  DTIME  on  the  order  of 
0.02  seconds  and  DTIME2  around  0.05  seconds  have  proved  to 
be  satisfactory.  A  TIMEF  of  0.5  to  0.7  seconds  is  satisfactory 
for  all  but  unusual  trajectories. 
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CARDS  9  and  10  contain  six  inputs  necessary  to 
start  any  trajectory  whose  initial  time  is  not  equal  to  zero. 
These  items  are  output  from  the  computer  program  and  are  input 


in  5 L 1 4 . 7  format 


VAR(l) 


VAR (2) 


VAR (3) 


VAR ( 4 ) 


VAR (5) 


VAR (6) 


is  the  x-location  of  the  store  center  of 
gravity  in  the  fudelage  coordinate  system, 
(feet) 

is  the  z-location  of  the  store  center  of 
gravity  in  the  fuselage  coordinate  system, 
(feet) 

is  the  store  pitch  angle  about  the  fuselage 
y-axis.  Positive  is  nose  up.  (degiaes) 
is  the  x-velocity  of  the  store  center  of 
gravity  in  the  fuselage  coordinate  system, 
(feet  per  second) 

is  the  z-velocity  of  the  store  center  of 
gravity  in  the  fuselage  coordinate  system, 
(feet  per  second) 

is  the  store  pitch  rate  about  the  fuselage 
y-axis.  Positive  is  nose  up.  (radians  per 
second) 
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Appendix  C 

Corputcr  P r o g r a n  Listing 

Explanation  o  f  Listed  Programs 

In  this  appendix  are  listed  the  main  program  and  all 
subroutines  with  the  exception  of  those  26  subroutines  taken 
from  Fernandes  (Ref  5).  Subroutine  ADAMS  was  taken  directly 
from  Goodwin  (Ref  6:76).  The  subroutine  FREAIR,  listed  here 
as  an  example,  was  composed  specifically  for  the  M-117  bomb 
and,  of  course,  will  be  different  for  other  stores.  Sub¬ 
routine  TBLNDC,  available  from  the  WPAFB  computer  library 
is  called  from  the  specific  FREAIR  listed  in  this  appendix, 
and  is  included  for  completeness. 

Changes  to  the  Fernandes  Subroutines 

Changes  were  necessary  to  the  Fernandes  subroutines  to 
pass  the  interference  coefficients  to  the  main  program  and 
to  avoid  extraneous  output. 

The  Fernandes  program  has  card  identifiers,  in  columns 
75  through  79,  sequenced  in  intervals  of  10.  Cards  listed 
below  whose  identifiers  end  in  the  number  "0"  are  to  replace 
those  in  the  Fernandes  subroutines  with  the  same  identifiers. 
All  other  cards  are  to  be  inserted  into  the  Fernandes  sub¬ 
routines  in  the  sequential  order  implied  by  their  identifiers. 
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